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^ : ABSTRACT 

^ ■ X-ray observations of the accreting X-ray pulsar 4U 1907+09, obtained 

^ ■ during February 1996 with the Proportional Counter Array on the Rossi X-ray 

• Timing Experiment (RXTE), have enabled the first measurement of the intrinsic 

; pulse period Pp^ise since 1984: Pp^ke = 440.34llg:E]}? s. 4U 1907+09 is in a 

' binary system with a blue supergiant. The orbital parameters were solved and 

■ - - ' this enabled the correction for orbital delay effects of a measurement of Ppuise 

obtained in 1990 with Ginga. Thus, three spin down rates could be extracted 
from four pulse periods obtained in 1983, 1984, 1990, and 1996. These are 
within 8% equal to a value of Ppuise =+0.225 s yr~^. This suggest that the 
pulsar is perhaps in a monotonous spin down mode since its discovery in 1983. 
Furthermore, the RXTE observations show transient ~18 s oscillations during a 
flare that lasted about 1 hour. The oscillations may be interpreted as Keplerian 
motion of an accretion disk near the magnetospheric radius. This, and the 
notion that the co-rotation radius is much larger than any conceivable value for 
the magnetospheric radius (because of the long spin period), renders it unlikely 
that this pulsar spins near equilibrium like is suspected for other slowing 
accreting X-ray pulsars. We suggest as an alternative that perhaps the frequent 
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occurrence of a retrograde transient accretion disk may be consistently slowing 
the pulsar down. Further observations of flares can provide more evidence of 
this. 

Subject headings: pulsars: individual (4U 1907+09) — stars: neutron — X-rays: 
stars — binaries 



1. Introduction 

4U 1907+09 is an X-ray pulsar powered by accretion of wind material from a blue 
supergiant companion star. It was discovered as an X-ray source by Giacconi et al. (1971) 
and has been studied using instruments onboard Ariel V (Marshall & Ricketts 1980, 
MR80), Tenma (Makishima et al. 1984, M84), EXOSAT (Cook & Page 1987, CP87), and 
Ginga (Makishima & Mihara 1992, Mihara 1995). MR80 first determined the orbital period 
of the binary at 8.38 days through an analysis of data taken over the course of 5 years 
(between 1974 and 1980) from a survey instrument on Ariel V with a net exposure time of 
about I year. A folded lightcurve of these data shows a pronounced primary flare and a 
dimmer and irregular secondary flare. M84 observed 4U 1907+09 with Tenma in 1983 and 
discovered the pulsar with a pulse period of 437.5 s. Also, through a good time coverage of 
the binary orbit, they were able to confirm the occurrence of two phase-locked flares. The 
secondary flare in the Tenma data appears as bright as the primary flare. CP87 discuss 
EXOSAT data with a small though reasonably uniform coverage of the orbit. They also find 
evidence for a phase-locked primary and secondary flare. Through combining Tenma and 
EXOSAT data, they were able to determine the binary orbit most accurately and found an 
eccentricity of O.lGlgji. A measurement of the pulse period revealed an average spin down 
rate of +0.23 s yr~^ since the Tenma measurement 270 d earlier. 4U 1907+09 was observed 
with Ginga in 1990. Makishima & Mihara (1992) report a cyclotron feature at 21 keV 
found during these observations. Mihara (1995) measured the pulse period using the Ginga 
data, without correcting for the binary orbit its value was determined to be 439.47 s. Sadeh 
& Livio (1982), using data from the HEAO A-1 instrument, report the occurrence of 15 ms 
oscillations during 4 out of 20 scans of the source each lasting about 10 s. The oscillation 
period was seen to change during each scan. 

In February 1996, 4U 1907+09 was observed with the narrow-field instruments on 
the Rossi X-ray Timing Explorer (RXTE). A previous paper (In 't Zand, Strohmayer & 
Baykal, 1997, ISB97) reported a frequent occurrence of complete dips in the pulsar signal 
with durations between a few minutes and l| hrs as found in these data. In the current 
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paper we report about the same data but concentrate on variability on time scales shorter 
than or equal to the pulse period. We discuss the pulse period and its history since 1983 
and present evidence for the occurrence of transient oscillations. These timing diagnostics 
enable one to address interesting questions about the presence of an accretion disk and 
the angular momentum transfer in the binary. Furthermore, we briefly discuss the general 
spectral trend as a function of pulse phase. 

2. Observations 

RXTE-PCA (Zhang et al. 1993) consists of five identical proportional counters that 
are co-aligned to the same point in the sky. Collimators limit the field of view (FOV) to 1 
square degree. The total geometric collecting area is approximately 6250 cm^, the effective 
sensitive photon energy range is from 2 to 60 keV. The satellite is in a low-earth orbit so 
that long observations may be broken up by earth occultations and passes through the 
South Atlantic Anomaly. 

RXTE-PCA observed 4U 1907+09 during 4 observation runs which are detailed in 
Table and total 79.4 ks. The mode of data collection for the four observation runs was 
event data with 64 energy channels resolution over the full bandpass and a time resolution 
of 16 fis. No detector identification bits were telemetred. At the time of the observations, 
the on-board gain correction algorithm was temporarily turned off (Jahoda et al. 1996). 
The gains of the 5 PCA detectors are slightly different. Whenever we quote intensities in 
a certain energy range, we refer to the energy range consistent with the channel setting of 
detector "PCUO" . We note that the gains of the other detectors are deviating by less than 
0.1 keV at the low end to at most 2 keV at the high end of the range. 

There are a number of known contributors to the X-ray background in our 
measurements: the cosmic diffuse background, the diffuse galactic ridge emission and the 
supernova remnant W49B. The latter two have not been accurately imaged yet in our 
bandpass. This forces us to employ a simplified background subtraction procedure which 
has been detailed in ISB97: the X-ray background is defined as the residual emission found 
during dip times when there is no apparent pulsed emission. Using this background implies 
two assumptions: if there is residual emission from 4U 1907+09 we disregard that and 
assume it is constant, and there is no other variable source in the background. The last 
assumption seems reasonable given the mentioned contributors, certainly on the time scale 
of the observations. With regards to residual emission from 4U 1907+09, we can only note 
that estimates of the different contributions are prone to systematic errors but suggest that 
residual emission is at a level less than 2% of the average 4U 1907+09 intensity in the 2 to 
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15 keV band (ISB97). 

3. The X-ray light curve 

Figure [l| shows the time history of the raw countrate of the RXTE-PCA observations 
on 4U 1907+09. Apart from the dips, the pulsar signal with a period of about 440 s is 
obvious. The pulse profile is very variable. It can change from pulse to pulse or even within 
one pulse. The time history shows a flare at Feb. 23.07 when the net 2 to 15 keV intensity 
rises up by an order of magnitude to about 0.1 Crab units. Given the newly determined 
orbital period (§ ^ this flare occurs 543.99±0.04 orbital periods after the secondary and 
544.44±0.04 orbital periods after the primary flare as observed with Tenma in August 
and September of 1983 (M84) which are the most recent timing reports of flares. The 
identiflcation of the flare on Feb. 23.07 with the secondary flare seems unambiguous and, 
therefore, we conflrm the recurrence of a phase-locked secondary flare. It is interesting to 
note that the orbital phase of the flare (see Table |l|) places it near to the time of apastron 
of the binary. The intensity of the flare is, in 2 to 30 keV, about twice as high as the one 
measured with Tenma (M84). This is not unexpected. The secondary flare has been seen 
to vary from insigniflcant to as bright as the primary flare (MR80). We note that the 
observed flare may actually be only part of a more extended secondary flare which could 
have been missed because of the sparsity of the RXTE coverage. The Tenma observation 
of the secondary flare suggests the flare may last 0.3 days, the associated uncertainty in 
the epoch has been included in the above-mentioned uncertainty of the number of orbital 
periods since the Tenma-observed flares. 

4. Timing analysis of the pulse signal 

An initial estimate of the pulse period was obtained by folding the 2 to 15 keV time 
history on a number of statistically independent trial periods (Leahy et al. 1983) in the 
range 430 to 450 s. Only data outside the intensity dips were used and photon arrival times 
were corrected into those for the solar system barycenter. The highest value was found 
for a period of 440.4 s (flgure Q). 

In order to accurately determine the pulse period as well as the binary orbit a set 
of 19 pulse arrival times was generated, one pulse arrival time for each XTE orbit when 
4U 1907+09 is not in a dip. This was done by folding the time history data into one 
average pulse for each RXTE orbit, folding all time history data into one master pulse, and 
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cross correlating the master pulse with each average pulse to find the pulse arrival time for 
each average pulse. Average profiles rather than individual pulses were used to minimize 
systematic effects due to the sometimes strong changes in the profile from pulse to pulse 
which are supposedly due to variability in the accretion rate (ISB97). 

As an alternative to control the pulse profile variability we have used the method 
of pulse wave filtering as proposed by Deeter & Boynton (1985, see also Boynton et al. 
1986). In this method pulse profiles are expressed in terms of harmonic series and cross 
correlated with the average pulse profile ('master pulse'). The maximum value of the cross 
correlation is analytically well-defined and does not depend on the phase binning of the 
pulses. The short term sharp fiuctuations of pulses are naturally filtered by a cut off of 
higher harmonics. The pulse arrival times obtained through this method gave statistically 
the same results as those obtained through the above-mentioned method. 

In order to increase the accuracy of the orbital solution (in particular the orbital 
period), the pulse arrival times were combined with the pulse delay-time data from Tenma 
observations as published graphically by M84. 

We modeled the data with an eccentric orbit with parameters that have not changed 
since the Tenma observations in September 1983 and with a pulse period which is constant 
throughout our RXTE observations. We determined the parameters by testing the model 
against the data using Pearson's test on a grid of parameter values sufficiently sampled 
to detect significant changes in and with a range to sufficiently enclose the allowed 
parameter values. The parameters are the pulse period Ppuise at the time of the RXTE 
observations, the orbital period Porb, the epoch of mean longitude 90 degrees (i.e., 
one quarter of the orbital period after the time of ascending node when the neutron star 
crosses the sky tangent plane through the barycenter moving away from the observer), 
the longitude of periastron uj (i.e., with respect to the ascending node), the eccentricity e, 
and the length of the projected semi- major axis ax sin i/c. The grid ranges and sample 
frequencies were found iteratively, going from a rough to a sufficiently fine grid. We choose 
this grid-search method to be able to determine the uncertainty of the solution with any 
arbitrary confidence level. The values were calculated in the pulse delay-time domain. 
On the one hand this enables the inclusion of the Tenma data that have not been published 
in the pulse arrival-time domain, on the other hand this precludes the determination of 
Ppuise during the Tenma observations. We fixed the error in the pulse delay times at 8 s. 
This value was suggested by the uncertainty of the harmonics in the pulse wave filtering 
analysis of the RXTE data and was the value M84 applied for the Tenma data. 

The results of modeling the timing data are presented in Table |^. The quoted error 
intervals are the projections of the 68% confidence level region onto each of the 6 parameter 
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axes. Since some parameters are very dependent, the actual solution is confined within 
a much smaller space than that of the 6-dimensional cube dictated by the quoted errors. 
The errors are, therefore, conservative. For the purpose of easy comparison with orbital 
solution by others, we have also included the single-parameter la errors in Table Figure ^ 
presents a plot of the orbital model versus the data in terms of delay times with respect to 
the binary barycenter. 

The orbital solution provides a means to correct the pulse period as determined with 
Ginga in 1990 (Mihara 1995) for the orbital motion of the neutron star in the binary system. 
The Ginga observations were performed for a duration of 14 hours (this is 0.07 in phase) 
at a mean time of MJD 48156.60. This is 0.10 orbital periods after the epoch of ascending 
node where the orbital Doppler shift is +0.28 s per pulse period with an uncertainty of 
about 0.01 s per pulse period. Therefore, the corrected pulse period is 439.19 s. 

We find a pulse period at the time of the RXTE observations of 440.34llo;oi7 s. This 
completes a total of four measurements of the pulse period and 3 measurements of the 
pulse period derivative since 1983. These values, tabulated in Table |^, imply an interesting 
finding: the pulsar appears to be spinning down with a close to constant derivative. The 
three values are within 8% of a mean value of Ppuise = +0.225 s yr^^. 

As mentioned, we assumed that over the course of 12 years there is no noticeable 
change in Porb, and the orbit inclination angle i. For Porb, we can confirm this with 
reasonable accuracy because the value we find is, within error margins, equal to what was 
found by CP87 between the Tenma and EXOSAT measurements (8.3745+0.0042 days). 
For uj, the uncertainty, typically 20 degrees, is simply too large to be at all able to measure 
likely values for the change (for a detailed analysis of apsidal advance in the similar system 
Vela X-1, see Deeter et al. 1987). The same argument holds for the sensitivity towards 
measuring changes in i. 

The photon arrival times were folded with the pulse period after correction for the 
earth's motion around the sun, RXTE's motion around the earth and the binary motion of 
the pulsar. The resultant pulse profile is presented in figure |^ for 6 photon energy bands 
up to 38 keV. Above this energy no pulsed emission was detected, as revealed by a Fourier 
analysis (figure ^). The pulse profiles are expressed in units of average 4U 1907+09 intensity 
per band. The background in each band was determined as discussed in ISB97 (see 
The pulse profile consists of two peaks with a deep and a shallow minimum in between that 
are close to 0.5 in pulse phase apart. The pulse profile looks similar to that observed more 
than a decade before with Tenma and EXOSAT in similar energy bands (M84 and CP87). 
It is rather insensitive to energy between 2 and 20 keV although subtle dependencies are 
noticeable, particularly in the second peak. There is a dramatic change of the pulse profile 
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at around 20 keV. Basically, the first pulse disappears above that energy and the shape of 
the other one changes. The modulation depth is roughly the same in all bands. The energy 
dependency of the pulse profile appears similar to that observed from Cen X-3 with Ginga 
(Nagase et al. 1992). 

5. Transient 18 s oscillations 

Visual inspection of the X-ray lightcurve on second timescales during the flare at Feb. 
23.07 suggests that there is a variety of non-Poissonian variability present during the flare. 
For example, figure |^ shows the lightcurve of a 500 second interval near the peak of the 
flare. Fluctuations in the countrate as large as several hundred counts s~^ are obvious in 
this lightcurve. To further quantify this short timescale variability (i.e., shorter than a 
pulse period) we computed Fourier power spectra for the time intervals encompassing the 
flare. The power spectrum of a 1024 second interval beginning 300 seconds prior to that 
shown in figure ^ is displayed in figure ^ There are a number of conspicuous peaks in 
this power spectrum in the range from 0.02 to 0.06 Hz. There is also clearly a broadband 
noise component increasing toward lower frequencies, as the mean power is significantly 
above 2 (the mean expected for purely poisson fluctuations) extending above 0.2 Hz. One 
of the most prominent peaks is at 0.055 Hz, or a period of 18.2 s. These 18 s oscillations 
can actually be seen with the eye in figure ^ To investigate further we performed an 
epoch folding period search in the vicinity of 18 s. The results are shown in figure ^, and 
there is an obvious peak centered at 18.2 s. To determine the average amplitude of the 
oscillations, we folded the 1024 s of data on the best-fit period of 18.2 seconds . We then fit 
a model including a constant countrate plus a sinusoid to the folded data. Figure ^ shows 
the resulting background subtracted lightcurve and the best-fit sinusoidal model. The 
per degree of freedom is formally a bit high, 1.8, but for estimating the amplitude of the 
oscillation the sinusoidal fit is sufficient. From this analysis we obtain an average oscillation 
amplitude of 4.4 ± 0.3%, where the amplitude is defined as the ratio of the sinusoidal to 
constant countrate components. We also investigated the dependence on photon energy of 
both the pulsation amplitude and pulse profile, but found no significant energy dependence. 

The 18 s oscillations are not persistent, rather, they are confined to an approximately 
1000 s interval centered near the peak of the flare. To investigate the transient nature of 
the 18 s oscillation we computed a dynamic power spectrum by calculating power spectra 
from 500 second intervals with a new interval beginning every 50 seconds. The resulting 
power spectra are not independent, since the data segments overlap, however, this method 
identifies the range of times for which the 18 s oscillations are present. The resulting 
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dynamic spectrum is shown in the top panel of figure |T^. The 0.055 Hz oscillations are 
clearly present from about 200-1200 seconds, and this is the only epoch in which we 
detected such oscillations. During the time that the oscillations are present there is no 
strong evidence for significant frequency drift, this gives a lower limit to the Q- value for the 
oscillation of Q = (1000s/18.2s) = 55.6 {Q is the so-called quality factor which measures 
the number of oscillations that will pass before a substantial fraction of the energy of the 
oscillator is dissipated, Q is inversely proportional to the bandwidth or spectral purity of 
the oscillation). 

To investigate the nature of the broadband noise we computed an average FFT power 
spectrum from two successive orbits during the fiare, for a total of 6250 seconds of data. 
The resulting average power spectrum is shown in figure |Tl]. We fit a power law model 
P = Ku^"' to 51 frequency bins in the 0.02 to 4 Hz frequency range. We did not fit 
below 0.02 Hz because there the power spectrum is dominated by the pulsed signal and its 
harmonics. The power law model with K = 0.41 ± 0.07 and a = 1.36 ± 0.06 provides an 
acceptable fit, with a x^/d.o.f. = 39/49 ~ 0.8. The integrated power from 0.02 to 4 Hz 
corresponds to an amplitude (rms) of about 8.2%. 

We note that we have searched for pulsations at higher frequencies and failed to find 
any. No pulsations are present like were reported by Sadeh & Livio (1982) at a period 
of around 15.3 ms. During the fiare we find an upper limit for the amplitude of 0.5% as 
compared to the amplitude of 12% at which pulsations were seen by Sadeh & Livio. 

6. Discussion 

Quasiperiodic oscillations (QPO) with frequencies in the 10 - 200 mHz range have 
been observed from 7 other X-ray pulsars. Among these, A0535-I-26 (Finger, Wilson, 
& Harmon 1996), X1626-67 (Shinoda et al. 1990), Cen X-3 (Takeshima et al. 1991), 
V0332+53 (Takeshima et al. 1994), X0115+63 (Soong & Swank 1989), and SMC X-1 
(Angelini 1989) all have QPO frequencies in the 0.062 to 0.1 Hz regime, similar to the 0.055 
Hz transient oscillation from 4U 1907-1-09 described above. The ~ 4% amplitude (rms) of 
the 18 s oscillation is also similar to the amplitudes of QPO from other X-ray pulsars (see 
Takeshima et al. 1994; and Angelini, Stella & Parmar 1989), however, QPO detected in 
these sources typically are broader, with Q-values of order a few compared to about 50 for 
the 18 s oscillation in 4U 1907-^09. 

To date, almost all models for QPO in X-ray pulsars postulate the existence of an 
accretion disk as the site of QPO production. In the case of QPO from A0535-I-26 this 
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is supported by the observed correlation of the QPO frequency and spin-up rate (Finger, 
Wilson, & Harmon 1996), strongly suggesting the transfer of angular momentum from 
accretion disk to neutron star. Thus, one possibility is that the 18 s oscillations reveal the 
presence of a transient accretion disk during the flare. 

Assuming that the 18 s oscillation in 4U 1907+09 is related to orbital motion, either 
via a beat frequency model (BFM, Lamb et al. 1985; Alpar & Shaham 1985), or a Kepler 
frequency model (KFM, van der Klis et al. 1987), then the inferred radius of material in a 
putative disk is = {GM / Ati'^u'^pq)^/^ ^ 1.2 x 10^ km, for a 1.4 Mq neutron star. This 
is nearly an order of magnitude smaller than the co-rotation radius -Rc = 9.6 x 10^ km for 
this object. We can rougly estimate the size of the neutron star magnetosphere from the 
expression for the Alfven radius assuming spherical accretion (e.g., Ghosh & Lamb 1991), 



Studies of the optical counterpart, discovered by Schwartz et al. (1972), reveal that the 
distance to 4U 1907+09 is between 2.4 and 5.9 kpc (Van Kerkwijk et al. 1989). Assuming 
the distance is 4 kpc, the average non-absorbed luminosity of 4U 1907+09 during non-dip 
and non-flare periods and within the studied energy range is 1 x 10^^ erg s~^ (ISB97). 
During the flare the luminosity increases up to 4 x 10^^ erg s~^. For a flare X-ray luminosity 
of 4 X 10^^ ergs/s, a magnetic field B — 2.b x 10^^ G and a stellar mass and radius of 
1.4 Mq and 10 km respectively, R^ k, 4300 km. Although this number is rather uncertain 
it suggests that the disk can certainly penetrate close enough to the star to account for 
the 18 s oscillations as an orbital phenomenon associated with an accretion disk. Although 
the accretion disk interpretation seems plausible, further observations, particularly in the 
vicinity of the flare will be required to confirm this interpretation. 

Long-term measurements of the behavior of the spin period for a selection of 15 
accreting X-ray pulsars (Bildsten et al., 1997) show that three known pulsars appear 
to exhibit systematic spin-down evolution for at least 5 years: GXl+4 {P/P ~90 yrs), 
4U 1626-67 (P/P 5 X 10^ yrs) and Vela X-1 {P/P Ri 6 x 10^ yrs). Generally, this behavior 
is attributed to low accretion rates that bring the magnetospheric radius R^a close to the 
co-rotation radius R^ (Ghosh & Lamb 1979a, 1979b; Wang 1987, 1995; Anzer & Borner 
1980; Lovelace, Romanova &: Bisnovatyi-Kogan 1995). However, for 4U 1907+09 this is an 
unlikely scenario because the magnetospheric radius is hard to bring out to such a large 
co-rotation radius unless either the magnetic field is of unlikely order 10^"^ G or the distance 



R. 



-m 



2.4 X 10^(M/1.4 Mq)^/^ X 

(P/2.5 X 10^'G)^/^ X 

(P/10^ cm)^*^/^ X 

(Lx/3.1 X lO^Vgs s-^)-^/^ km. 



(1) 
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is of order 0.5 kpc which is at least 5 times lower than studies of the optical counterpart 
suggest. 

The putative disk supports the notion that the magnetospheric radius is substantially 
smaller than the co-rotation radius. Moreover, a disk can provide an alternative 
explanation to the observed spin down if this transient disk is rotating in an opposite 
sense to the pulsar spin. Such a disk can provide the necessary torque to explain the 
spin down rate. If all accreted mass supplies its angular momentum at a radius the 
expected torque on the neutron star is Nchar = r]M(GMR(iy^'^ where < 1 is the duty 
cycle of the transient disk. If we assume that all the liberated potential energy of the 
accreted mass during the flare is transformed into radiation and the flare luminosity is 
3.2 X lO^^Dlp^ erg s~\ then A^^hax = 2.5 x lO'^^riDl^^R]/^ . For R^ = 1.2 x 10^ km this 
becomes A^char = 8.6 x lO^^TyD^pc g cm^. The observed absolute value of the torque is 
A^obs = 27r /|z>| where I is the moment of inertia and z> = —3.7 x 10~^^ Hz s~^ (equivalent 
to P = +0.225 s yr~^). If 7 = 10^^ g cm^ (generic value for a neutron star of radius 10 km 
and mass 1.4 M0), iVobs = 2.3 X 1032 g cm2. If iVobs = iVchar and 2.4 < Dkpc < 5.9, then 
0.008 < 1] < 0.05. Therefore, the putative retrograde transient accretion disk can provide 
the torque to spin down the pulsar if it lasts on average on the order of a few percent of the 
time. We observed a duty cycle of about 1000 s out of an orbital period of 8.4 d which is 
an order of magnitude too small. However, the coverage of our observations is limited and 
also the flare might, averaged over many orbits, have a larger duty cycle. This is confirmed 
by Tenma observations which suggest a duration of 0.3 d for one flare (M84). We conclude 
that the putative retrograde accretion disk could possibly supply enough negative torque to 
spin down the pulsar. 

The suggestion of a retrograde accretion disk being responsible for an extensive 
spin down period has recently been revisited in the case of the disk-fed AXP GXl+4 by 
Chakrabarti et al. (1997) in order to have an elegant explanation for a positive correlation 

between luminosity and negative torque. Om need for a retrograde accretion disk in the 
wind-fed 4U 1907+09 is motivated by a co-rotation radius which is clearly much larger than 
the magnetospheric radius. 

4U 1907+09 is in several ways similar to Vela X-1: the spin periods are both of order 
several hundred seconds (440 s for 4U 1907+09 and 283 s for Vela X-1), the orbital periods 
are of order 8 days (8.38 d for 4U 1907+09 and 8.96 d for Vela X-1), and both have recently 
been found to show long periods of systematic spin down (e.g., Bildsten et al. 1997 for 
Vela X-1). Previously, the occasional spin-down periods in Vela X-1 have been attributed to 
the presence of a disk with an inner edge close to the co-rotation radius (Nagase 1989), with 
the note that this implies relatively large magnetic fields (surface fields of order 10^^ G). 
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We suggest that, since this reasoning goes awry in the case of 4U 1907+09, this conclusion 
may be invahd for Vela X-1 as well. An important difference between both systems is that 
Vela X-1 shows much more spectral variability and pronounced iron lines. This may be due 
to a difference in angle of sight of the binary orbit, a difference in the spatial distribution of 
the wind from the companion star, and/or a somewhat different eccentricity. 

There is, as compared with some other AXPs, moderate spectral variability within 
the pulse profile. Basically the pulse profile shows two peaks below and one peak above 
~20 keV. This is somewhat similar to what has been observed from the much more 
luminous disk-fed AXP Cen X-3 and may be associated with two magnetic poles which have 
either different physical circumstances or are viewed on at different angles by the observer. 
A modehng of the pulse profile, similar to studies performed by Buhk et al. (1995) and 
beyond the scope of this paper, may provide more definite conclusions about this. 

Since its discovery by MR80, the occurrence of two flares per orbit has been difficult 
to reconcile with the identification of the companion as a supergiant high- mass star (e.g.. 
Van Kerkwijk et al. 1989, ISB97). It has always been thought that these stars do not have 
a circumstellar disk like is presumed for Be stars. However, recently evidence is mounting 
that in hot supergiants in the upper part of the H-R diagram axial symmetry may play an 
important role (e.g., Zickgraf et al. 1996 and references therein). Also, recently another 
wind-fed AXP with a supergiant companion has been shown to consistently display two 
flares per orbit: GX 301-2 (Pravdo et al. 1995; Chichkov et al. 1995; Koh et al. 1997). 
This again supports the notion that the wind from the companion is not isotropic but that 
the mass flux is enhanced along the equatorial plane. The two flares then may be caused 
by the neutron star traversing this plane twice per orbit in a sufficiently inclined orbit. It is 
interesting to note in this context that GX 301-2 (Sato et al. 1986) and 4U 1907-^09 have 
the highest eccentricity of the group of supergiant high-mass XRBs, perhaps these systems 
are not yet in a tidal equilibrium. Perhaps even the spin orientation of the neutron star 
differs by more than 90° from the orbital motion orientation. However, at the moment this 
is mere speculation. 



7. Conclusion 

We have determined the pulse period of 4U 1907-1-09 at a recent epoch and find 
that the pulsar has spun down on average by 0.225 s yr~^ since the pulsar's discovery 
in 1983. Three measurements of the spin down rates during intermediate time interval 
spanning between 1 year and 6 years are consistent within 8% of this value. This suggests a 
remarkable monotonous spin down trend during about 12 years. Furthermore, we find the 
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occurrence of 18 s oscillations for 10^ s of a flare with a high Q value hinting at an accretion 
disk. The indication that the magnetospheric radius is much smaller than the co-rotation 
radius, the observed long and constant spin down trend and the occurrence of transient 
oscillations suggest that possibly a recurrent transient accretion disk counter rotates the 
neutron star and, through transfer of angular momentum, slows the pulsar down. This 
suggestion needs to be confirmed by at least detailed multiple observations of flares which 
are almost certainly recurring every orbital period. 
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Table 1: RXTE observation log of 4U 1907+09 



Observation run 


1 


2 


3 


4 


Dates (Feb. 1996, U.T.) 


17.73-18.15 


19.67-20.10 


21.60-22.08 


23.02-23.42 


Orbital phases^ 


0.55-0.60 


0.78-0.84 


0.01-0.07 


0.18-0.23 


Orbital phases'^ 


0.85-0.90 


0.08-0.14 


0.31-0.37 


0.48-0.53 


Exposure time (s) 


19442 


19911 


20261 


19775 


Time span (s) 


36093 


37601 


41658 


33457 



"Based on epoch of maximum distance between neutron star and earth 
''Based on epoch of periastron 



Tabic 2: The binary orbit and pulse period of 4U1907-I-09 from 1983 Tenma and 1996 
RXTE-PCA measurements 



Parameter 


Symbol 


Value 


68% 


Single 








confidence 


parameter 








region 


1(7 error 


Orbital period 


-forb 


8.3753 days 


+0.0003 
-0.0002 


0.0001 


Eccentricity 


e 


0.28 


+0.10 
-0.14 


0.04 


Orbital epoch 




MJD 50134.76 


+0.16 
-0.20 


0.06 


Longitude of periastron 




330 degrees 


+20 
-20 


7 


Projected semi-major axis length 


Ox sin ijc 


83 h-s 


+4 

-4 


2 


Pulse period 


-fpulse 


440.341 s 


+0.012 
-0.017 


0.006 



Tabic 3: History of Ppuisc iiicasureiiieiits for lU 1907+09 



Mean time 




SateUite 


Reference 


Ppulse 


Derivative^ 


Date 


(MJD) 






(s) 


(s yr-i) 


Aug./Sept. 1983 


45576 


Tenma 


M84 


437.483+0.004 




May/June 1984 


45850 


EXOSAT 


CP87 


437.649+0.019 


0.22+0.03 


Sep. 1990 


48156.6 


Ginga 


Mihara 1995 


439.19+0.02*^ 


0.244+0.005 


Feb. 1996 


50134 


RXTE 


this paper 


440.341+0.014 


0.212+0.004 



"The pulse period derivative is calculated from the difference in Ppuise with respect to the previous 

observation in this table. 

''This value was corrected for delays from binary motion in the present work. The uncertainty is an estimate. 
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Fig. 1.— Time histories of the raw RXTE-PCA 2 to 15 keV intensity of 4U 1907+09, the 
bin time is 32 s. The orbital phases are with respect to the epoch of the maximum distance 
between the pulsar and the solar barycenter 
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Fig. 2. — Periodogram of the 2 to 15 keV time history data after correction to solar system 
barycenter but before correction to binary system barycenter 
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Fig. 3. — The delay time versus time for the RXTE measurements (+ symbols) and for the 
Tenma measurements folded into the RXTE time of observation (x symbols). The phases 
indicated on top are with respect to the time of maximum distance between the pulsar and 
the solar barycenter. The solid line shows the model for the eccentric orbit specified in 
Table |^. The vertical dashed lines refer to the inferred times of periastron. All data points 
are repeated modulo Porb- The rms of the residuals is 6.6 s 
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2- 3 keV 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
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Fig. 4. — The background-subtracted lightcurve folded with the pulse period for 6 bandpasscs 
between 2 and 38 keV. The bandpasses are indicated in each panel as well as the statistical 
1(7 error in the upper left corner (from top to bottom these errors are equal to 0.045, 0.007, 
0.009, 0.010, 0.018 and 0.038). The phase offset is arbitrary. The unit of intensity is the 
average intensity per bandpass. Only data for 'quiet' periods is used when the source is 
neither dipping nor flaring above 500 c s~^. The net exposure time for these lightcurves is 
29.6 ks out of a total of 79 ks 
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Fig. 5. — Periodogram of the > 40 keV time history data after correction to solar barycenter 
but before correction to binary barycenter 



-21 - 




Fig. 6. — An expanded view of the observed X-ray time profile during the secondary flare in 
the full PCA bandpass. Time is measured from Feb. 23, 1996, 01:56:33 UTC 
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Fig. 7. — The power spectrum of a 1024 s piece of the hghtcurve starting 300 s before the 
start time indicated in figure | 
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Fig. 8. — Periodogram near 18 s for the flare data 
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Fig. 9. — Folded background-subtracted lightcurve, folded with a period of 18.2 s, and a 
sinusoidal fit to the data 
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Fig. 10. — The top panel shows a contour map of the dynamic power spectrum. The bottom 
panel shows the accompanying lightcurve. The vertical dashed lines indicate the time frame 
of figure H 
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